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Abstract. This article summarises the status of the solar neutrino oscillation phenomenology at 
the end of 2002 in the light of the SNO and KamLAND results. We first present the allowed areas 
obtained from global solar analysis and demonstrate the preference of the solar data towards the 
Large-Mixing-Angle (LMA) MSW solution. A clear confirmation in favor of the LMA solution 
comes from the KamLAND reactor neutrino data. The KamLAND spectral data in conjunction with 
the global solar data further narrows down the allowed LMA region and splits it into two allowed 
zones - a low Am^ region (low-LMA) and high Am^ region (high-LMA). We demonstrate through a 
projected analysis that with an exposure of 3 kton-year (kTy) KamLAND can remove this ambiguity. 
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1. The Neutrinos from the sun 

Solar neutrinos are produced via the reaction 

4p i/e + 2e+ + 2iye + 28MeV (1) 

The above process occurs through two main cycles of nuclear reactions - the pp chain 
(CNO cycle) which is responsible for 98.5% (1.5%) of the energy. There are eight dif- 
ferent types of neutrino fluxes, named according to the parent nuclei of the decay chain 
which generates it. The pp chain gives rise to the neutrinos pp, pep, hep, ''Be, 
while the neutrinos ^^N, ^'^0,^^F are generated through nuclear reactions forming the 
CNO cycle. The solar neutrino fluxes are calculated by the so called "Standard Solar Mod- 
els" (SSM) among which the most extensively used are the ones due to Bahcall and his 
collaborators [1]. The flux predictions from the SSM are robust. Different solar models 
agree to a very high degree of accuracy ( to within 10%) when the same input values of 
the parameters are used and also demonstrate striking consistency with helioseismological 
measurements. The pp neutrinos are mainly responsible for solar luminosity and the SSM 
prediction for the pp flux is least uncertain. The prediction for the ^B neutrino flux is most 
uncertain stemming from the uncertainties associated with the cross-section of the reaction 
Be{pj)^B producing these neutrinos. 
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2. Solar Neutrino Experiments 



The pioneering experiment for the detection of solar neutrinos is the ^"^Cl experiment in 
Homestake which started operation in 1968 ^. It utihses the reaction [4] 

i^e +^'^ CI ^^"^ Ar + e- . (2) 

The threshold for this is 0.814 MeV and hence it is sensitive to the and ^Be neutrinos. 

Three experiments SAGE in Russia and GALLEX and its updated version GNO in Gran- 
Sasso Underground laboratory in Italy uses the reaction [5] 

Ue + ^^Ga ^^Ge + e" (3) 

for detecting the solar neutrinos. This reaction has a low threshold of 0.233 MeV and the 
detectors are sensitive to the basic pp neutrinos. 

The radio chemical experiments ^"^Cl and ^^Ga experiments are sensitive to only Uf. and 
can provide the total solar i^e flux. 

The first real time measurement of the solar neutrino flux was done by the Kamiokande 
imaging water Cerenkov detector, located in the Kamioka mine in Japan [6]. It was sub- 
sequently upgraded to SuperKamiokande - a same type of detector but with much larger 
volume increasing the statistics [7]. The neutrinos interact with the electrons in the water 
via 

e~ ^ e~ + Vx (4) 

This reaction is sensitive to all the three neutrino flavours. However the and i^r react 
via the neutral current which is suppressed by a factor of 1/6 compared to the Ve inter- 
action which can be mediated by both charged and neutral currents. The recoil electron 
energy threshold in Kamiokande was 7.5 MeV which could be reduced to 5 MeV in Su- 
perKamiokande. Thus both the detectors are sensitive mainly to the neutrinos. 

The Sudbury Neutrino Observatory (SNO) experiment also uses a Cerenkov detector but 
containing heavy water {D2O). The deuterium in heavy water makes it possible to observe 
solar neutrinos in three different reaction channels [8,9] 

Ue + d^p + p + e- (CC) 
I'x + e' ^ i'x + e~ (ES) 

Vx+d^n+p + Vx (NC) 

The charged current (CC) reaction is exclusive for v^, The electron scattering (ES) reaction 
is same as in SK. The unique feature of SNO is the neutral current (NC) reaction which 
is sensitive to all the three flavours with equal strength. For both CC and ES reactions 
the final state electrons are directly detected through the Cerenkov hght emitted by them 
which hits the PMTs and an event is recorded. For the NC reactions the final state neutron 



^For recent reviews on solar neutrino experiments see [2,3]. 



2 



can be captured (i) by another deuteron (ii)by capture on CI in an NaCl enriched heavy 
water (iii)by ^He proportional counters. For both (i) and (ii) the nuclei after capturing the 
neutrons emits single and multiple gamma rays respectively which compton scatters the 
electrons in the medium. The Cerenkov light produced by these electrons will produce 
an event. Therefore if the NC events are due to (i) and (ii) above then they cannot be 
disentangled from the CC and ES events Exclusive detection of the neutrons produced in 
the NC event is possible for the process (iii). 



2. 1 The total solar neutrino flux 



The ratio of the observed solar neutrino rates to the SSM predictions are presented below. 



experiment 



obsvd 



CI 0.337 ± 0.029 
Ga 0.553 ± 0.034 
SK 0.465 ± 0.014 
SNO(CC) 0.349 ±0.021 
SNO(ES) 0.473 ±0.074 
SNO(NC) 1.008 ±0.123 



The declared SNO NC data is due to neutron capture on deuteron. Therefore the CC,ES 
and NC events cannot be separated on an event by event basis. For extracting the separate 
rates from the entangled data sample one one needs to assume an undistorted flux as an 
input. The SNO rates quoted in the above Table are obtained under this assumption [9]. In 
all the above experiments observed v,, flux is less than the theoretical predictions implying 
disappearance of the solar UgS. On the other hand for the SNO NC data the observed rate 
agrees to the theoretical prediction. Since the NC is sensitive to and Vr as well this 
indicates that the are reappearing as v^s and/or z/t-s. 



2.2 Information on Direction and Energy 

Apart from providing a measurement for the total solar neutrino flux the real time measure- 
ments can also provide information on direction and energy of the incoming neutrinos. The 
electron scattering reaction used in SK and SNO has excellent directional sensitivity. In 
fact through this reaction the Kamiokande experiment first demonstrated the solar origin of 
the neutrinos. The left panel in figure 1 plots the number of events observed in SK against 
the cosine of the angle with the sun's direction. There is a clear peaking towards cos 9q 
= 1.0 [10]. The statistical capacity of SK allows it to make an energy wise binning of the 
data and present the recoil electron energy spectrum. The right panel in figure 1 shows 
the data/SSM as a function of the recoil electron energy with a 5 MeV threshold. The plot 
exhibits a flat recoil electron energy spectrum consistent with no spectral distortion [10]. 

The panels 1 and 3 in Figure 2 show the number of CC+ES+NC+background events 
in SNO as a function of direction and energy respectively. The solid lines in the figures 
show the Montecarlo simulated CC,ES and NC events. The ES events show the strong 
directional correlation with sun as in SK. The CC events has an angular correlation 1 — 
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Figure 1. The left panel shows the number of events observed in SK as a function of 
cos ^0 while the right panel shows the energy spectrum of the recoil electrons observed 
by SK. 



0. 340056*0 while the NC events have no directional dependence as the produced gamma 
rays do not carry any information of the incident neutrino. The electron produced in the CC 
reaction being the only light particle produced in the final state has a strong correlation with 
the incident neutrino energy and thus can provide a good measurement of the energy 
spectrum. The recoil electron spectrum from the ES reaction is softer. The NC events due 
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Figure 2. The SNO spectrum vs cos ^0, radial coordinate and recoil electron energy. 
Also shown is the MonteCarlo simulated events. 



to capture on deuteron has a peaked distribution around the energy of 6.25 MeV. These 
probability density functions are used to perform a maximum likelihood fit to the data. 
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2.3 Variations with time 



SK already has enough statistics to divide their data into both energy and zenith angle bins. 
The latest SK data has been presented with a binning of eighth energy bins and each energy 
bin contains data subdivided into seven zenith angle bins. [11]. This binning enables one to 
measure the day time and night time fluxes separately and study for any possible day/night 
asymmetry. Since in day time the neutrinos do not pass through the earths matter and in 
night time they traverse the earths matter an observed day /night asymmetry would be clear 
indicator of earths matter effect modifying the neutrino fluxes. The day/night asymmetry 
measured in SK and SNO are stiU not at a statisticaUy significant level. 



2.4 Evolution of the solar neutrino problem 

Before the declaration of the SNO data there were two aspects of the solar neutrino prob- 
lem. The first one is - in aU the experiments the observed flux was less than SSM pre- 
diction. The second one was the 'problem of the ^Be neutrinos'. We note that among the 
pre-SNO experiments SK and its predecessor Kamiokande is sensitive to the neutrinos. 
The ^''Cl experiment is sensitive mainly to the ^B and ''Be neutrinos while the sensitivity 
of the ^^Ga experiments amount to pp, ^B and ^ Be neutrinos. Combining the ^B flux 
measured in SK with the CI data resulted in no room for the ''Be neutrinos. Similarly 
the expected pp flux in Ga consistent with solar luminosity plus the ^B flux observed in 
SK indicated a negative flux of ^Be neutrinos in Ga. The solar physics could not explain 
this preferential vanishing of ''Be flux over ^B flux as in the pp chain ''Be comes before 
^B and any mechanism that reduces the ^ Be flux would eventually reduce the ^B flux. 
Therefore the answer was sought in the properties of neutrinos and neutrino flavour con- 
version was considered as the most promising candidate for the solution. SNO provided 
the compelling evidence. 

The SNO CC reaction is sensitive only to while the ES reaction in both SK and SNO 
is sensitive to both Ue and v^Ivt. Therefore a higher ES flux as compared to the CC flux 
wiU imply the presence of 1/^/1/,- in the solar flux. Combining the SK ES and SNO CC 
results one gets 

<fES - (l>cxf = 0-57 ± 0-17 X lO'^/cmVsec (5) 

This is a 3.3 fT signal for Vf, transition to an active flavour (or against transition to solely 
a sterile state). 

The NC reaction is sensitive to all the three flavours with equal strength resulting in a 
greater sensitivity to the neutral current component in the solar flux. Comparing the NC 
and CC data from SNO one gets 

- (l^cxP = 3-41 ± 0.65 x lOVcmVsec (6) 

In two circumstances we can have the CC and NC rates equal to each other. Either when 

there is no flavour conversion or for flavour conversion to a purely sterile state which 
does not interact with the detector. The observed CC/NC difference rules out both these 
possibilities at 5.3cr. 
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3. Two Flavour Oscillation 



If neutrinos have mass then the flavour eigenstates Ve and y^/i'T are different from the 
mass eigenstates 1^1,^2 and related as 



ye \ _ ( COS 9 sm.9\ f vi 
J \ — sin ^ cos 6 I \ V2 



(7) 



where 6 is the mixing angle in vacuum. This leads to neutrino oscillation in vacuum [12]. 
Then the survival probabihty that a Vf, remains Vf. after traveling distance L in vacuum is 

P^^^^ = 1 - sin226'sin2(1.27Am2L/£;) (8) 

Am^ = — m\. The term containing Am^ is the oscillatory term resulting from coher- 
ent propagation of the mass eigenstates. 

In matter, only z/g's undergoes Charged current interaction giving rise to an matter in- 
duced mass term of the form \/2GF-^e- This changes the mixing angles as 

Am^sin2^ 

tan26'M = ^ 7- (9) 

Am2cos2(9-2\/2GFnei^ 

is the electron density of the medium and E is the neutrino energy. Eq. (9) demonstrates 
the resonant behavior of Om- Assuming Am^ > the mixing angle in matter is maximal 
(irrespective of the value of mixing angle in vacuum) for an electron density satisfying, 

2V2GF Ue.resE — Am^ cos29 (10) 

This is the Mikheyev-Smirnov-Wolfenstein (MSW) effect of matter-enhanced resonant fla- 
vor conversion [13]. 

The most general expression for i^,, survival probability in an unified formalism over the 
mass range 10~^^ — lO^'' eV^ and for the mixing angle 9 in the range [0,7r/2] is [14] 

Pee = P0P© + (l-P0)(l-Pffi) 



+2v'P0(l - P0)P®(1 - P©) cos^ (11) 

where Pq denotes the probability of conversion of i^^ to one of the mass eigenstates in the 
sun and P© gives the conversion probability of the mass eigenstate back to the i^e state in 
the earth. All the phases involved in the Sun, vacuum and inside Earth are included in ^. 
Depending on the value of Am^ /E one has the following three limits 

(i) in the regime Am^/E ~ 5 x 10^^*^ eV-^/MeV matter effects inside the Sun suppress 
flavor transitions while the effect of the phase ^ remains. This is the vacuum oscillation 
limit. 

(ii) For Am?/E ~ 10~^ eV^/MeV, the total oscillation phase becomes very large and 
the cos^ term in Eq. (11) averages out to zero signifying incoherent propagation of the 
neutrino mass eigenstates. This is the MSW limit. 

(iii) For 5 X 10-1" eV^/MeV ~ Am^/E ~ lO'^ eV^/MeV, both matter effects inside the 
Sun and coherent oscillation effects in the vacuum become important. This is the quasi 
vacuum oscillation (QVO) regime. 
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Figure 3. The allowed regions in Am^ -tan^ 9 plane from data on total solar rates in 
the left panel and from data on total solar rates and SK spectrum data in right panel. 



Next we present the results obtained by performing a x^-analysis of the solar neutrino 
data. The procedure followed can be found in [15-17]. In figure 3 we present the al- 
lowed regions obtained from analysis including the total fluxes measured in CI, Ga, SK 
and SNO [17]. There are basically five regions which are allowed - the small mixing angle 
region (SMA), the large mixing angle high Arn^ regions (LMA), the large mixing angle- 
low Am^ regions (LOW), the vacuum oscillation regions symmetric about tan^ 6 = 1.0 
and the quasi vacuum region between the LOW and the vacuum oscillation regions. 

In the scond panel of figure 3 we present the allowed areas after including the SK spec- 
trum data with the total rates data [17]. The SMA region and large part of the vacuum 
oscillation region are seen to have been washed away with the inclusion of the SK spec- 
trum data. 

In the left panel of fig. 4 we show the dependence of the probabilities on energy. In the 
SMA and the VO oscillation regions the probability has a non-monotonic dependence on 
energy whereas in the LMA and LOW regions the survival probability does not have any 
appreciable dependence on energy beyond 5 MeV which is the threshold for SK and SNO. 
Thus these regions are favoured by the flat SK spectrum. 

In the right panel of figure 4 ( from [ 1 8]) we present the allowed regions obtained from an 
global analysis [18-21] of the total rates observed in CI and Ga (SAGE, GALLEX and 
GNO combined rate), the 1496 day SK zenith angle energy spectrum data and the recent 
SNO data of combined CC+ES+NC+background data in 17 day and 17 night energy bin. 
After the inclusion of the SNO results 

(i) LMA emerges as the favoured solution. 

(ii) The LOW region now appears only at 3(t. 

(iii) Values of Am^ above 3 x 10^^ eV^ are seen to be disfavored. 
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E (in McV) tan 9 

Figure 4. The left panel shows the dependence of the probability on energy. The right 
panel gives the allowed area after including the latest SNO spectrum data including the 
NC events. 



(iv) The QVO and VO solution are not allowed at 3(t. 

(v) Maximal mixing (6 — 7r/4) is disfavored at 3.4ct. 

(vi) SMA solution is disfavored at 3.7a 

(vii) The Dark Side {9 > tt/A) solutions disappear. 

Apart from including the SNO spectral data the figure 4 also has the new 1496 day SK 
zenith angle spectrum data. The data reveal a flat zenith angle spectrum. The inclusion of 
this data rule out the higher part of the LOW solution beyond 1.5 x 10~^ eV^ for which 
peaks in the zenith-angle spectrum are expected [22]. 

4. KamLAND 

KamLAND is a 1000 ton liquid scintillator neutrino detector situated at the former site 
of Kamiokande [23]. It looks for iJ^ oscillation coming from 16 reactors at distances 81 
- 824 km. Most powerful reactors are at a distance ~ 180 km. The detection process is 
i7f,p e+n. The e+ produced annihilates with the electrons in the medium to produce 
prompt photons. The neutrons get absorbed by the protons in the medium to produce de- 
layed photons. The correlation of time,position and energy between these two constitute 
a grossly background free signal. The survival probability relevant in KamLAND is the 
vacuum survival probability (cf. eq. 8) summed over the distances from all the reactors. 
The average energy and length scales for KamLAND are ^ 3 MeV, L ~ 1.8 x 10^ m 
which makes it sensitive to Am^ ^ 1.6 x 10^^ eV^ which is in the LMA region. The 
figure 5 shows the probabilities for KamLAND for the average distance of 180 km and 
solar neutrinos for Am^ (=6 x 10^^ eV^) and tan^ 9 (=0.5) [24]. Whereas for the solar 
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probabilities the Ato^ dependence is completely averaged out in KamLAND the probabil- 
ity exhibits a L/E dependence which gives it an unprecedented sensitivity to determine 
Am^ . 




E (MeV) 

Figure 5. The left figure shows the variation of the survival probability with energy 
for KamLAND and for the solar neutrinos. The right figure gives the allowed region 
from an analysis of KamLAND rate and global solar data. The allowed region from 
only solar analysis is shown also by the outer dashed line. 

The solar best-fit predicts a rate in KamLAND O.GS^Q gg (3a) [25] while the observed 
rate is 0.611 ± 0.094 corresponding to 145 days of data. Thus KamLAND confirms the 
LMA solution. In figure 5 we show the allowed area obtained from a x^-analysis of Kam- 
LAND rates data and global solar data [26]. The contour obtained from only solar analysis 
is also drawn (dashed lines) for understanding the role played by KamLAND . The inclu- 
sion of the KamLAND rates data gives a lower bound of Am^ > 4 x 10~^ eV^. The other 
parts of the parameter space allowed from the solar analysis are not constrained much. 

Figure 6 shows the allowed areas obtained from a analysis of the only Kam- 
LAND spectrum [23]. For an average energy of 5 MeV and a distance 180 km a Am? of 
7 X 10" 5 eV^ corresponds to the oscillation wavelength (A) ^ the distance traveled(L) 
and Pee ~ 1, for Am^ of 3.5 x lO"'"^ eV^ \ = 2L and Pee = 0.0. Again for Am^ = 
1.4 x 10~^ eV^ A = L/2 and Pee — 1. Since the KamLAND spectral data corresponds to 
a peak around 5 MeV islands around the first and the third Am^ are allowed whereas the 
middle An? is disfavoured as is seen from the panel 1 of figure 6. 

In the second panel of figure 6 we show the allowed area from KamLAND spectrum and 
global solar data obtained through a analysis [26]. Inclusion of the KamLAND spectral 
data splits the allowed region into two zones at 99% C.L. low-LMA (LMAl) and high- 
LMA (LMA2). LMA2 has less statistical significance (by w 2a) The global best-fit comes 
at Arn^ = 7.17 X 10"^ eV^ and ta,i? 9 = 0.44 in low-LMA. LOW region is disfavoured at 
5a. Maximal mixing although allowed by the KamLAND spectrum data gets disfavoured 
at 3.4ct by the overall analysis. 

In the left panel of figure 7 we explore through a projected analysis of 1 kton year 
simulated spectrum the potential of KamLAND in distinguishing between the two allowed 
areas. 

The allowed Am? ranges around both low-LMA and high-LMA zones decrease in size. 
Since the solar data favours the low-LMA zone the allowed areas become more precise for 
spectrum simulated at the low-LMA zone while ambiguity between the two zones remains 
for high-LMA spectrum. In the right panel of 7 we show that a higher statistics (3 kton 
year) from KamLAND can remove this ambiguity almost completely determining Am^ to 
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Figure 6. The left panel shows the allowed regions in the Am^ — tan^ 6 plane from 
analysis of KamLAND spectrum data. The right panel shows the allowed region from 
a global analysis of KamLAND spectrum and solar neutrino data. 



within 6%. 

Analysis of solar and KamLAND data has been done in the realistic three neutrino sce- 
nario in [27]. In this case the non-zero value of the mixing angle 613 ^ connecting the 
atmospheric sector and the solar sector modifies the allowed area in Am? — tan^ 9 pa- 
rameter space. A third solution at a Am^ higher than the high-LMA zone gets marginally 
allowed for the three flavour case. As 613 is increased the high-LMA zones tend to disap- 
pear. 



5. Conclusions 



The solar neutrino research began with a motivation of understanding the sun through the 
neutrino channel. Over the years it metamorphosed into a tool of unraveling the fundamen- 
tal properties of the neutrinos. At the end of 2002 the status of the solar neutrino oscillation 
phenomenology can be summarised as 

• Comparison of SNO CC and SNO NC signifies neutrino flavour conversion at 5.3(T. 

• Rules out transitions to pure sterile states at 5.3(T. 

• LMA is the favoured solution to the solar neutrino problem with 3 x 10^^ eV^ < Am^ < 
3 X 10"'* eV2 and 0.25 < tan^ 9 < 0.87 at 3ct. 

• KamLAND confirms LMA. 

• Best-fit Ato^ shifts from 6 x 10^^ eV^ as obtained from global solar analysis to 
7.2 X 10-5 eV2. 

• There is no significant change in best-fit 6 (tan^ 9 = 0.4). 
•LMA region splits into two parts at 99% C.L.. 

• 3ct allowed range after the first KamLAND data is 4.96 x lO'^ eV^ < Am^ < 2 x lO""* 
eV2 and 0.27 < tan^ 9 < 0.88. 



^this is at present limited by the CHOOZ data to sin^ 613 < 0.04. 
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Figure 7. The 1 kTy (left panel) and 3 kTy (right panel) projected contours of Kam- 
LAND and solar data. The points at which the spectra are simulated are shown by bold 
dots. 



•Transitions to a mixed state is still allowed with < 13%(52%) sterile mixture at la{3a) 
[24]. 

At this juncture the emerging goals in solar neutrino research are 

(i) precise determination of the oscillation parameters; 

(ii) to observe the low energy end of the solar neutrino spectrum consisting of the pp,CNO 
and the "^Be line and do a full solar neutrino spectroscopy. 

As far as precision determination of Am^ is concerned KamLAND has remarkable sensi- 
tivity and 3 kton year KamLAND data can remove the ambiguity between the two presently 
allowed zones completely. Even before that a day-night asymmetry > 1% in SNO can ex- 
clude LMA2 [28]. However the constraining power of KamLAND for 9 is not as good 
being limited by the 6% systematic error and also by the fact that in the statistically sig- 
nificant regions the observed KamLAND spectrum corresponds to a peak in the survival 
probability where the 9 sensitivity is very low. At the present value of the best-fit param- 
eters in the low-LMA region a new KamLAND like reactor experiment with a baseline of 
^ 70 km will be sensitive to the the minimum in the survival probability increasing the 9 
sensitivity by a large amount [29]. The upcoming Borexino experiment should see a rate 
0.64 ± 0.03 (Icr) and no day-night asymmetry. It cannot differentiate between the two 
LMA regions. However it can provide a measurement of the ^Se flux coming from the 
sun. Real time measurement of the pp neutrino flux is the target of the so called "LowNU" 
experiments like XMASS, HELLAZ, HERON. CLEAN, MUNU, GENIUS using ly - e 
scattering and LENS, MOON and SIREN using charged current reactions [30]. All these 
experiments involve development of new and challenging experimental concepts and re- 
search work is in progress for evolving these techniques. 
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